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Confining ultracold gases in cavities creates a paradigm of quantum trapping potentials. We 
show that this allows to bridge models with global collective and short-range interactions as novel 
quantum phases possess properties of both. Some phases appear solely due to quantum light-matter 
correlations. Due to global, but spatially structured, interaction, the competition between quantum 
matter and light waves leads to multimode structures even in single-mode cavities, including delo¬ 
calized dimers of matter-field coherences (bonds), beyond density orders as supersolids and density 
waves. 


Ultracold atoms trapped in optical lattices (OLs) en¬ 
able to study quantum many-body phases with unde¬ 
niable precision and target problems from several disci¬ 
plines [T]. Such optical potentials can be complicated, 
but are prescribed, i.e., they are created by external 
lasers and are not sensitive to atomic phases. This lim¬ 
its the range of obtainable states. Self-consistent light- 
matter states can be obtained, when scattered light mod¬ 
ifies the trapping potential itself. This was achieved by 
trapping a Bose-Einstein condensate (BEG) inside an op¬ 
tical cavity HH, which dramatically enhances the light- 
matter coupling, thus making the influence of reemission 
light comparable to that of external lasers. Such “dynam¬ 
ical potentials” [5] enabled the structural Dicke phase 
transition and a state with supersolid properties [2]. A 
key effect observed so far is the dynamical dependence 
of light intensity (potential depth) on the atomic den¬ 
sity. Although, the light becomes dynamical, its quan¬ 
tum properties are still not totally exploited as works 
on atomic motion in quantum light were limited to few 
atoms [6H9]. Effects in dynamical potentials are anal¬ 
ogous to semiclassical optics, where atomic excitations 
are quantum, while light is still classical. As the light 
and BEG are quantum objects, the quantum fluctuations 
of both were studied [Mill], however, the fundamental 
reason of the structural phase transition can be traced 
back to the dynamical self-organization predicted [12] 
and observed m with thermal atoms and classical light. 
Eor single-mode cavities, dynamical light-matter cou¬ 
pling was shown to lead to several effects [HHISI yet 
to be observed with bosons and with non-interacting 
fermions [TM2T] . Recently an optical lattice inside a cav¬ 
ity has been realized [^ . Multimode cavities extend the 
range of quantum phases further [S] . 

We show that, even in a single-mode cavity, the quan¬ 
tum potential [elETlEHl leads to significant many-body 
effects beyond semiclassical ones. Multimode spatial pat¬ 
terns of matter fields arise due to symmetry breaking 
resulting from the competition between imposed global 
light structure and standard local processes (tunneling 
and on-site interactions). We demonstrate that the effi¬ 
cient competition is achieved due to the ability to struc¬ 
ture the global interaction at a microscopic scale con¬ 



FIG. 1: Cold atoms trapped in an optical lattice subject to a 
quantum potential created by the light inside a single-mode 
cavity. The unsharp potential contour schematically depicts 
quantum fluctuations of light, which induce the light-matter 
correlations. The cavity can be a standing- or traveling-wave. 
Different colors represent atoms corresponding to different 
light-induced spatially structured atomic modes. 

sistent with the lattice period. Such a competition in 
turn leads to novel many-body states, which are not lim¬ 
ited to density-induced orders as in previous studies, but 
also represent long-range patterns of matter-field coher¬ 
ences (bonds [in]), leading, e.g., to far delocalized dimers, 
trimers, etc. Importantly, we prove that our approach 
bridges models with global collective and short-range in¬ 
teractions, as new quantum phases possess properties 
of both, going beyond Dicke, Lipkin-Meshkov-Glick [33] 
and other simple spin-1/2 models. Recently, nontrivial 
spatial patterns were obtained with classical atoms and 
light [3D]- Our work will assist to extend such efforts 
in the interdisciplinary field of optomechanics towards 
quantum multimode systems m- The mechanisms we 
suggest, provide a general framework and a new set of 
tools, inaccessible in setups using classical OLs. This will 
strongly expand applications in quantum simulations. It 
will allow exploring fundamental issues concerning emer¬ 
gence of multimode generalizations of strongly corre¬ 
lated phases, such as gapped superfluids [33] and density 
waves [33] as well as their interplay, giving rise to quan¬ 
tum solids [33] ■ The light-induced structure is similar to 
multi-component nonlinear sigma models ubiquitous in 
analog models of high-energy [SD] [36], condensed mat¬ 
ter pIllSH], and relativistic [39] physics. Dimer phases 
can be used as building blocks for quantum spin-liquids 
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simulations [lU] . 

We consider atoms trapped in an OL inside single¬ 
mode cavity with the mode frequency lOc and decay rate 
K in off-resonant scattering (see Fig. [^. The pump light 
with the amplitude flp (in units of the Rabi frequency) 
and frequency Wp (Ac = Wp — Wc) illuminates atoms in 
a plane transverse to the cavity axis, but not necessarily 
at 90°. The atoms couple with cavity mode via the ef¬ 
fective coupling strength g 2 = g^p/(2Aa), where g is the 
light-matter coupling coefficient and is the detuning 
between the light and atomic resonance ^51 ITTl 02] ■ This 
can be described by the Hamiltonian H = H'’ +'H°‘ , 

where Ti,^ is the regular Bose-Hubbard (BH) Hamilto¬ 
nian [4^47j . The light is described by 7^“ = hiOcO^a and 
the light-atom interaction is |28) : 

= g*2aF^ + g^b) F ( 1 ) 

with F = D + B. D = Jjjfij is the diagonal coupling 

of light to on-site densities, B = j) + h.c.) is 

the off-diagonal coupling to the inter-site densities reflect¬ 
ing matter-field interference, or bonds mi. The sums go 
over illuminated sites Ng, {i,j) denotes nearest neighbour 
pairs. The operators (a) create (annihilate) photons 
in the cavity, while (bi) correspond to bosonic atoms 
at site i. 77“^ is a consequence of the quantum potential 
seen by atoms on top of BH model given by a classical OL 
with the hopping amplitude to and on-site interaction U. 

The spatial structure of light gives a natural basis to 
define the atomic modes, as the coupling coefficients Jij 
can periodically repeat in space and are calculated from 
the Wannier functions [13] see m- The symmetries bro¬ 
ken in the system are inherited from such a periodicity: 
all atoms equally coupled to light belong to the same 
mode, while the ones coupled differently belong to differ¬ 
ent modes. We define operators corresponding to modes 
^ , where 

Dtp — Jr)'p ; with Np — ^ ^ Tii , (^) 

with V = I] {b\bj+h.c.).{i) 

Thus, we replaced the representation of atomic operator 
T as a sum of microscopic on-site and inter-site contri¬ 
butions by the smaller sum of macroscopically occupied 
global modes with number density, and bond, 
operators.The structures of density and bond modes can 
be nearly independent from each other. To be precise, 
for the homogeneous scattering in a diffraction maxi¬ 
mum, Jij = Jb and Jjj = Jd, one spatial mode is 
formed. When light is scattered in the main diffraction 
minimum (at 90° to the cavity axis), the pattern of light- 
induced modes alternates sign as in the staggered field, 
Jij = Jjp = (—1)-^ Jb and Jjj = {—lyjD- This gives 
two spatial density modes (odd and even sites) and, as 


we will show, four bond modes. The density and bond 
modes can be decoupled by choosing angles such that 
Jd = 0 (by shifting the probe with respect to classical 
lattice thus concentrating light between the sites and as¬ 
suring the zero overlap between Wannier and mode func¬ 
tions) or Jb =0 [H]. Beyond this, additional modes get 
imprinted by pumping light at different angles such that 
each i?-th site or bond scatters light with equal phases 
and amplitudes. This generates multimode structures 
of R density modes [321133] and 27? bond modes. The 
prominent example of self-organization [HIIMII] is a 
special case of two density modes, while macroscopic ef¬ 
fects related to the higher density modes and any bond 
modes have not been addressed so far. 

In general, the light and matter are entangled [6l|27l 
IMH53] . In the steady state of light, it can be adiabatically 
eliminated and the full light-matter state can be then 
reconstructed as we show in m- The effective atomic 
Hamiltonian [ilIil[2H] is 

nls=n^+ ^{F^F + FFy, (4) 

where (?eff = ^c| 52 p/(Ac-l-K^). A key physical processes 
is that the ground state is reached [i.e. the energy @ 
is minimized], when the system adapts (self-organize) in 
such a way that the light scattering term is maximized 
for ^eff < 0, and minimized for g^g > 0. New terms be¬ 
yond BH Hamiltonian give the effective long-range light- 
induced interaction between density and bond modes: 


f^f+ff^ = 




(5) 


= symmetry 


where 

broken by the light modes imprints the structure on the 
interaction of atomic modes. 

Fundamentally, Eq. (101 displays the link between 


global interactions and the interaction resembling typi¬ 
cal short-range one (usually appearing between the sites, 
while here between the modes). Thus, the resulting quan¬ 
tum phase will have properties of both collective and 
short-range systems. In this language, it is the term 
iV^dd-^even f^at is responsible for supersolid properties 
of the self-organized state (The standard supersolidity 
appears due to interaction.) Our general ap¬ 

proach enables to go far beyond typical global models 
(e.g. Dicke and Lipkin-Meshkov-Glick [31]) due to spa¬ 
tial structuring the global interaction thus assuring its 
effective competition with the short-range ones, even in 
a single mode cavity in contrast to [S] I35H35] . The bond 
interaction can be easily identified as pseudo-spin inter¬ 
action via the Schwinger mapping m- In addition, some 
components can be non-Abelian. 


We decompose (10) in mean-field contributions and 
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fluctuations: 

+ = (f1')F + (F)F^ + ^fI'F. (6) 

The last term SF^F originates from the quantum light- 
matter correlations, underlying the quantumness of OL. 
Other terms originate from the dynamical but classical 
light, when the semiclassical approximation oF^ = {a)F^ 
holds. Decorrelating operators at different sites, we ob¬ 
tain a mean-field theory that has nonlocal coupling be¬ 
tween the matter modes and is local in fluctuations. For 
SD^D these reduce to on-site number fluctuations. Im¬ 
portantly, this corresponds to the purely light-induced ef¬ 
fective on-site interaction of atoms beyond the standard 
BH term. For SWB, light-matter correlations include 
radically new terms beyond BH model: fluctuations of 
the order parameter and density coupling between neigh¬ 
boring sites, which appear due to two and four point 
quantum atomic correlations m- In contrast to previ¬ 
ous works, we will show non-negligible effects due to such 
terms, putting forward the quantumness of OL. 

When the ground state of is achieved by maximiz¬ 
ing scattering (^eff < 0), a strong classical light emerges 
and small fluctuations can be neglected. In principle, 
even in the strong-light case, the light quantumness can 
play a role, because self-organized states can be in a su¬ 
perposition of several patterns and different light ampli¬ 
tudes are correlated to them [B]. Nevertheless, in a real¬ 
istic case with dissipation, the system quickly collapses 
to one of the semiclassical states [54] . We will show that 
quantum fluctuations play a key role in the opposite case, 
where scattering is minimized [g^g > 0). Here no classi¬ 
cal light builds up and light fluctuations design the emer¬ 
gence of novel phases. To underline key phenomena, we 
will consider cases with either density or bond modes. 

Scattering at 90°, one breaks the translational symme¬ 
try. Hence, the system can support density waves (DWs) 
and novel bond orders. The simultaneous occurrence of 
SF and DW orders is a supersolid (SS) phase [34]. SS and 
DW have been predicted due to classical maximized scat¬ 
tering [Min]- We show that for weak pump, DWs and 
SSs with only small density imbalance appear at half¬ 
integer Ailing, together with usual MI and SF EUllT]. 
In contrast, above the threshold > U/2, we find 

that DWs and SSs with maximal imbalance are favored, 
while usual MI and SF are completely suppressed [^ITT] . 

For minimized scattering (geff > 0), the classical light 
cannot build up at all, and quantum fluctuations take 
the leading role [Fig[^a)]. At fixed density per site, the 
quantum light-matter fluctuations effectively renormal¬ 
ize the on-site interaction from U to U + 2gegJ^ in each 
component. Thus, changing the light-matter coupling, 
one can shift the SF-MI transition point. This occurs 
because the light-induced atomic fluctuations now enter 
the effective Hamiltonian, and these fluctuations need to 
be suppressed to minimize the energy. This favors MI 
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FIG. 2: Light scattering from on-site densities, (a) Mini¬ 
mized scattering highlighting the quantumness of OL. Total 
order parameter Ei/; = |'!^+P -f I'l/'-P, where 'ip± are the SF 
order parameters of two modes, white lines correspond to 
Mis; dashed lines are boundaries of SF with only two non- 
negligible Fock components. The system is homogeneous, 
p+ = p- and ipg = ij}-. White points are the MI-SF tran¬ 
sition points without cavity light, the SF-MI can be signifi¬ 
cantly shifted by the quantum OL. (b) Strongly interacting 
phase diagram for multiple number of modes R (created for 
different pump angles) at half-filling. Quantum phases have a 
period of density-density correlations q (in units of lattice pe¬ 
riod), DWq,! (DW 5 ,s), l{s) denotes large (small) density im¬ 
balance DW, A is the SF fraction [47l[55l[56j. DW and SF or¬ 
der depend on the effective light-matter interaction strength 
(?eff and R. Horizontal lines denote phase boundaries be¬ 
tween quantum phases. Parameters: (a) scattering at 90°, 
(7eff = IQU/Ns, the boundaries are for geft = IQU/Ns and 0; 
Jd = 1.0, Jb — 0, Ns = 100, z = 6. (b) scattering at vari¬ 
ous angles defining the mode number R, Jb = 1.0, Jb = 0, 
to =0, Ns = 100. 


state for U smaller than that without cavity light, extend¬ 
ing MI regions [Figj^a)]. For incommensurate fillings, 
SF survives but with smaller (suppressed) fluctuations 
as well: For convenience, we also plot the boundaries 
where the superfluid ground state is composed of mainly 
the two lowest occupation Fock states depending on the 
filling factor (i.e. components with higher occupations 
are negligible). Note that with cavity light, the state be¬ 
comes gapped with respect to adding more than one ex¬ 
citation, thus minimizing fluctuations. However there is 
no phase transition to this peculiar superfluid state [47j . 
Moreover, in a quantum OL, atoms can potentially en¬ 
ter MI even without atomic interaction. This provides 
absolute control on DW order formation. In analogous 
fermionic systems DWs are relevant for the stability of 
superconducting phases Ezj. 


Scattering at angles different from 90° creates more 
than R = 2 atomic modes (42 [43], and the light-imposed 
coefficients are Jjj = JDX(.j)i where for traveling waves 
x(j) = exp (j27rj/F). Now, multiple terms N^N^, in 
Eq. (10) become important. We present a phase dia¬ 
gram for various scattering angles (inducing R modes) for 
strong on-site interaction with maximal light scattering 
[Figj^b)]. The i?-mode induced pattern competes with 
on-site interaction modifying the density distribution. 
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Therefore, multiple DWs of period R can coexist with 
SF, forming multicomponent SSs. Surprisingly, at half¬ 
filling for i? > 2 odd, SS exists. As |(?eff| changes from 
zero, checkerboard insulators form for even R while dif¬ 
ferent kinds of SS exist for odd R. The on-site interaction 
limits atomic fluctuations producing gapped SF compo¬ 
nents when SS exists. As R increases, additional DWs 
with different periods and amplitudes emerge. These 
form unstable mixed state configurations (T) for i? = 3 
and stable multicomponent SS for R = 5. These occur 
in between regular SS phases with DWs of period R and 
different amplitude. As the light induces atomic fluc¬ 
tuations, these generate competition between small and 
large amplitude DWs until saturation in the SF compo¬ 
nent occurs. The system reaches a configuration, similar 
to the maximal imbalance DW state described above for 
scattering at 90°. 

We find a novel phase transition, when light scatter¬ 
ing from the bonds at 90° {Jjj+i = (—1)-^Js, Jd = 0) 
is maximized (^eff < 0). Even in the absence of on-site 
interaction, a transition from normal SF to the super¬ 
fluid dimer (SFD) state appears. SFD is a SF state in 
which the complex order parameter has alternating (zero 
and non-zero) phase difference between pairs of sites, and 
its amplitude is modulated as well. This occurs because 
of the competition between the kinetic energy BH terms, 
which promote a homogeneous SF, with the light-induced 
interaction that favors SF components with alternating 
phases across every other site [Figj^a)]. The phase of 
light interference flips from bond to bond [such that the 
phase difference at neighboring bonds is tt, FigUa)]. In 
the Hamiltonian, this corresponds to alternating signs 
in front of matter-field coherences between the neigh¬ 
boring sites (i.e. products of complex order parameters, 
in mean-field treatment). To minimize the energy (and 
maximize the light scattering) the quantum matter fields 
self-organize such that the matter-field phase difference 
between neighboring sites flips as well to compensate 
for the imposed light-held phase flip. The dimer con- 
hgurations have high degeneracy as the full many-body 
ground-state is composed of several equivalent arrange¬ 
ments of the phase pattern in space [47] . The phase di¬ 
agram is shown in Figj^b). Moreover, in the presence 
of on-site interaction the system supports a transition to 
the supersolid dimer (SSD) state with modulated densi¬ 
ties [Figj^c)] . The on-site interaction suppresses atomic 
huctuations, while as light scattering gets optimized, the 
density is unable to lock in a homogeneous density pat¬ 
tern leading to additional density imbalance. Therefore, 
the phase modulation and density modulations coexist si¬ 
multaneously while atoms retain mobility preventing the 
stabilization of an insulating phase. Note, that mul¬ 
timode bond structures can have very nontrivial spatial 
overlap and dimers (and their multimode generalizations 
as trimers, tetramers, etc., which can be obtained for 
tilted pump angle) extend over many sites demonstrat- 
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FIG. 3: Emergent dimer phases and quantum-light-induced 
supersolids due to scattering from bonds, (a) Dimer struc¬ 
ture for maximized scattering: matter-field coherence com¬ 
pensates for imposed light-held patterns, (b) Phase diagram 
for the phase difference between dimers when light scat¬ 
tering is maximized without on-site interaction, (c) Phase 
diagram for the difference in order parameters ipA/B of the 
dimers Aip = — |V’sPl/2 when light scattering is max¬ 

imized with on-site interaction, (d) Density wave order pa¬ 
rameter Ap = Ipa —Pfl|/2 for minimized scattering. The den¬ 
sity components Pa/b correspond to atomic populations of ef¬ 
fective light-induced modes. Regions with Ap ^ 0 correspond 
to supersolid phases; Ap = 0 corresponds to superfluid. Pa¬ 
rameters: (b) U = 0, (c) peff = —25U/Ns, (d) peflf = 25U/Ns; 

Jd = 0, Jb = 0.1, Ns = 100, and 2 = 6. 

ing the interplay of global and short-range properties. 
Dimer states can be used as fundamental units to en¬ 
gineer Hamiltonians of quantum spin-liquid states [40] . 
As W B processes enter directly in the effective Hamilto¬ 
nian with couplings depending on the light geometry and 
pump amplitude, this makes feasible to achieve analogous 
physics relaxing the constraint on very low temperatures 
and large on-site interactions based on second order ex¬ 
pansion effective Hamiltonians [T]. 

Importantly, we prove that there is a SS to SF tran¬ 
sition that is solely driven by quantum correlations for 
minimized light scattering (^eff > 0),Figj^d). This 
occurs because the terms due to light-matter correla¬ 
tions in B are not shadowed by semiclassical ef¬ 
fects, as there is no classical light build up. Two-point 
tunneling correlations introduce new terms in the BH 
model m, which couple densities at neighboring sites 
only: J2{i, fiihj, producing a DW instability even with¬ 
out strong light. Density imbalance is energetically 
favoured and the atoms condense in a nearest-neighbour 
density pattern, while additional terms in favor 

atomic quantum fluctuations competing with the on-site 
interaction. Thus, short range processes induced by the 
quantumness of light induce the transition. Such a direct 
density coupling corresponds stronger to the typical su¬ 
persolidity scenario, which is under active search [331 EH] ■ 

In experiments with homogeneous BECs both k and 
the long-range interaction rate are of the same order (ei- 
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ther MHz or kHz [HI]), while the depletion rate [SI] 
of the ground state can be made smaller by choosing 
the detuning such that |k/Ac| <1. In our lattice case, 
the effective light-matter coupling coefficient g rescales 
as gJD,B, and similarly to the homogeneous systems the 
evolution can be faster than depletion. 

In conclusion, we have shown that quantum optical 
lattices offer a new tool to engineer nonlocal many-body 
interactions with light-induced structures. These interac¬ 
tions can break symmetries by design and imprint a pat¬ 
tern that governs the origin of many-body phases. This 
effectively bridges physics of long-range and short-range 
interactions. The light and matter are entangled, forming 
non-trivial light-matter correlated states. We suggested 
how to generate not only multimode density patterns, but 
nonlocal patterns of matter-filed coherences as well (in 
particular, delocalized superfluid and supersolid dimers). 
Some of the states appear solely due to quantum fluctu¬ 
ations of light and matter, where no classical light can 
build up. A pathway to realize our proposal is to combine 
several recent experimental breakthroughs: a BEC was 
trapped in a cavity, but without a lattice EHl], and de¬ 
tection of light scattered from ultracold atoms in OL was 
performed, but without a cavity [HU IH3] • Very recently 
an optical lattice in a cavity became a reality [22] . Based 
on off-resonant scattering and thus being non-sensitive 
to a detailed atomic level structure, our approach can be 
extended to other arrays of natural or artificial quantum 
objects: spins, fermions, molecules (including biological 
ones) [HI] , ions [HH] , atoms in multiple cavities [HH] , semi¬ 
conductor [Hz] or superconducting qubits [68] . 

The work was supported by the EPSRC 
(EP/I004394/1). 
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Classical optical lattice model. The Bose-Hubbard 
Hamiltonian [151 US] is 


+ h.c) - - 1), 

(ij> i i 

(7) 


where {i,j) refers to nearest neighbour pairs, (hi) cor¬ 
respond to creation (annihilation) operators of bosonic 
atoms at the site i and the atom number operators are 
hi = The tunneling amplitude of the bosons is to, 

the on-site interaction is U and the chemical potential is 
/i. The effective parameters of the Bose-Hubbard Hamil¬ 
tonian with the cavity field can be calculated from the 


Wannier functions and are given by 

to = y w(x-Xi)(V^ - VbL(x))w(x-Xj)d”x, 

( 8 ) 

Jij = j w(x — Xj)u*(a;)itp(x)r(;(x — Xj)d"x, (9) 

where Uc^p{x) are the cavity and pump mode functions 
and u>(x) are the Wannier functions. The classical opti¬ 
cal lattice potential is given by Vol(p^)- Typical values 
for the amplitudes of couplings for the standing-wave po¬ 
tential Vol{x) = Vo sin^(27rx/A) with Vo « bEn are: 
to ~ O.IEr, \Jg^\ ~ 0.07,|J^| « 0.85, where Er is the 
recoil energy. This has been calculated using real Wan¬ 
nier functions finding the maximally localised generalised 
Wannier states of the classical optical lattice using the 
MLGWS code [15] . 

Light-induced interaction decomposition. New 

terms beyond BH Hamiltonian give the effective long- 
range light-induced interaction between density and bond 
modes: 

F^F + FF^ = E 

+ ,tpJB ,ip' + <^-C-)S,pS^, 

+ {JD,ip'^B,Lpi h-C-c)N^pS^i 

+ {Jb,pJd,ip' C.c)N^pS^i. ( 10 ) 

We separate light matter-correlations and dynamical 
terms in performing on-site mean-field. The D 
term can be expanded as 


b^b + bb^ « Y,{j*d,pJd,p' + c.c.){n^.){2n^- 

- (E))> 



+ sb'^b 

( 11 ) 

sb^b = 2J2\Jd,^\"snI 

ip 

( 12 ) 

i€ip 

( 13 ) 


where {N,p) = Pi i® mean number of atoms in 

the mode ip and pt = {hi) is the mean atom number at 
site i. The first term in Eq. 0 is due to the dynamical 
properties of the light field, these terms exhibit nonlo¬ 
cal coupling between light-induced modes. The terms 
in (121 signifies the light-matter correlations and contain 
the effect due to quantum fluctuations. 

The WB terms can be expanded as: 


B^B + BB^ « Y.{rB,^JB,^'+c.c.){S^,){2S,p-{S^)) 

-V SB'^B, (14) 

5B^B = 2Y,\Jb,v?5SI, (15) 
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(ij>6V 

+ E (^I^I(^)' + + ‘^n.blb^ + blh^ 

{i,j,k)Gv 

- (Pkh +^IK + + h.c.)), (16) 


where {i,j,k) refers to i,j nearest neighbours and k is 
a nearest neighbour to the pair {i,j). The first term in 
(141 is due to the dynamical properties of the light field 
and (16) are due to light-matter correlations. These are 
basically all the possible 4 point correlations and tun¬ 
neling processes between nearest neighbours, as higher 
order tunneling processes have much smaller amplitudes. 
In principle, additional terms might be considered in the 
above expansion and their generalization is straight for¬ 
ward by removing the restriction over the sums beyond 
nearest neighbours. The on-mode light matter correla¬ 
tions + £>1,13^ reduce to on-site covariances per 

mode. In the above we have decorrelated products of op¬ 
erators at different sites such that + {vj)£ — 


where ^ or i) are combinations of 6^ or b opera¬ 
tors at a given site. The expectation value of the bond op¬ 
erators reduces to (S^) = +ip]ipi), which 

is the sum of products of order parameters at nearest 
neighbour sites in the light-induced mode tp. In princi¬ 
ple additional incident or output light-modes might be 
considered by trivially adding a sub-index and sum over 
additional incident or output modes, this is equivalent 
to include additional light induced bands in the effective 
Hamiltonian. The technique can be fully transferred to 
multimode cavities and the coupling of additional cavities 
in the system. 


Effective Hamiltonians. The effective Hamiltonian 
considering only diagonal coupling {D'^ D + DD^) is 


^eff — ~^o E(^1^7 + 
{iJ) 

+ E [E 


(17) 


1) 2^gff|T£) (pi Piflr^ 




/i(p — P Pefi'kjDjtp: 

= U + 2geff\JD^ip\'^, 


(18) 

(19) 


with rjD,^ = + c.c.)(7Vp/) - |Jd,(pP 

and CD,(p = + c.c.){N^>){N^)/2 - 

\Jd,ip\'^ Pi - The many-body interaction U^p and the 
chemical potential fx^p inherit the pattern induced by the 
quantum potential that depends on light-induced mode 
structure given by ip. 

In the case of only off-diagonal bond scattering {B^B + 


BB^) , we have 

^eff — ^ ^ tipSip -\- I ,CyC I 

( 20 ) 

+ - phA , (21) 

i ^ ' 

t(p — tg gei^PBjcpt AA 

with 77 b,( p = + c.c.)(V) and cb,p = 

J2p>iJB,ipJB,ip' + c.c.)(,S,p')(,S,p)/2. The effective tunnel¬ 
ing amplitude tp couples the SF components of all the 
light-induced modes p. The terms due to <55'^ induce 
nontrivial coupling between nearest neighbour sites and 
lead to the formation of a density wave instability with 
more than one light-induced mode. 

The effective coupling parameters riB/B,ip depend im¬ 
plicitly on all the light induced modes giving nonlocal 
coupling via the expectation values of the operators 
and Np. 

Homogeneous scattering at diffraction marxi- 

mum. In the specific case of a single light-induced 
mode component in the diffraction maxima, we have: 
77 b = 2J'jj[NsP — 1/2) and tjb = ‘2zJgNs\ip\'^ with 
p = {fii) and Ip = (bi) for all sites and z is the coor¬ 
dination number. This gives the effective tunneling am¬ 
plitude tp = to — 2 z( 7 eff JgA^s|7/p, the effective chemical 
potential /ip = p — — Ap— 1) [where we have 

added all onsite density terms], and the effective inter¬ 
action strength Ucp = U + 2geffj‘l,. Density-dependent, 
but classical light was previously shown to strongly mod¬ 
ify the standard phase diagram of Mott insulator (MI) 
- superfluid (SF) transition [TU [IS], if plotted via the 
chemical potential p. If the phase diagram is plotted 
via the density, rather than chemical potential, most of 
such modifications due to the classical light become invis¬ 
ible. Therefore, in this paper we chose to plot the phase 
diagrams via the density to focus more on the effects 
produced by the quantum light. For example, at fixed 
density per site, it is the quantum light-matter correla¬ 
tions that effectively renormalize the on-site interaction 
from U to U + 2gesJl) [see Fig. |^a)]. Thus, changing the 
light-matter coupling, one can shift the SF-MI transition 
point dne to essentially the quantumness of light. 

Choosing geometry lu, one can suppress the density 
scattering (Jb = 0) and have all Jb’s equal. As no 
symmetry is broken, the bond self-organization does not 
emerge, but another semiclassical effect arises: tunnel¬ 
ing is enhanced (suppressed) for maximum (minimum) 
scattering. This modifies the phase diagram (Fig. |^), 
because of nonlinear coupling of the SF order parameter 
ip = {bi) to the tunneling amplitude to, which renormal¬ 
izes to to - 2zgesJlNs\'ip\‘^. 

A Peculiar Superfluid: the quantum superpo¬ 
sition (QS) state. In the presence of cavity light, the 
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FIG. 4: Modifications of quantum phases due to quan¬ 
tum and semiclassical effects for homogenous scat¬ 
tering. (a) Phase diagram in terms of SF order parame¬ 
ter ‘ip at fixed density for density-induced scattering. For 
minimal (maximal) scattering, MI boundaries (white lines) 
become extended (shortened) with respect to the transition 
point without cavity light (white point). This corresponds to 
suppression (enhancement) of quantum-light-induced atomic 
fluctuations. The behaviour of boundaries of the gapped QS 
state is similar (black dotted lines). The white dotted line 
is for gapless QS without cavity light, (b) Phase diagram 
for off-diagonal bond-induced scattering. The processes and 
lines are similar to (a), but arise due to renormalization of the 
tunnelling amplitude resulting from semiclassical light scat¬ 
tering. Parameters: (a) geft = 25U/Na, the boundaries are 
for get! = 25U/Ns, 0 , and —12.5U/Nb’, Jd = 1 . 0 , Jb = 0 ; (b) 
Qes = l.QU/Na, the boundaries are for gefi = l.GU/Ns, 0 , and 
— l.QU/Na', Jd = 0, Jb = 0.05; (a,b) — 100 and z = 6 

(3D). 


FIG. 5: Emergence of density wave and supersolid 
phases in the diffraction minima of light with den¬ 
sity coupling for maximized light scattering from the 
atoms (^eff < 0). (a) Phase diagram for the total order pa¬ 
rameter of two density modes T,ip = \ip+\'^ + IV’- 1^) where ip± 
are the SF parameters of two modes. The DW order param¬ 
eter is Ap = |p+ — p-l/2 when Ap 7 I 0 there is DW order 
( p± are the mean atom numbers per site in each mode). 
Dashed lines at half-integer fillings correspond to DW insu¬ 
lators (EV> = 0, Ap 7 I 0), SS state occurs, when Titp ^ 0 and 
Ap 7 I 0. White solid lines at integer fillings correspond to 
MI; dotted lines are the boundaries between SS and SF. (b) 
Phase diagram where only giant DW and SS exist above the 
quantum critical point for light-matter coupling (no SF or MI 
phases appear). The white points correspond to the classical 
optical lattice SF to MI transition.The boundaries of MI, SF, 
and MI can be tuned by the pump-cavity detuning. Parame¬ 
ters: (a) geff = -0.5t//iVs, (b) fleff = -1.25[//iVs, = 100 , 

z = 6 (3D). 


gap opens in the SF state because the effective chem¬ 
ical potential depends on the density and renormalizes 
to pL — geff</|)((2A^s — l)p — !)• Without tunneling for 
geS > 0 , the energy required to add a particle on top of 
the ground state is AEqs{p) = Up +geffJ‘^{2Nsp+l) for 
incommensurate fillings between MI regions with fillings 
n and n -I-1, with 


- ^ fH- \ 

^ 2g,sJlNa \U V 


(23) 


and on-site number fluctuations A(fii) = (p—n)(l—p+n). 
For MI regions at commensurate density p = n, the gap is 
AEuiin) = Un + gesJf){2Nen + l) and A{ni) = 0. This 
means that only occupations of the lowest particle-hole 
excitations are allowed in between MI lobes. There¬ 
fore in mean-field approximation this peculiar SF state 
is made of the quantum superposition of only two lowest 
Fock components that would satisfy the constraint on the 
density, instead of a coherent superposition of all possible 
fillings. Therefore, we use a term quantum superposition 
(QS) state for such a peculiar SF with a gap. The usual 
MI physics with renormalized interaction strength occurs 
at commensurate fillings. Similarly to SF-MI transition 
point, the boundaries of QS can be tuned (stabilized) 
due to the quantum correlations as U is renormalized 
[FigQa)]. The quantum superposition state can be un¬ 
derstood as the SF state of minimal atomic fluctuations. 
Note that a gap opens in the spectrum of excitations in 


the system, but because the particle filling is not com¬ 
mensurate with the lattice the system can’t form an insu¬ 
lator and the best alternative energetically is the QS. As 
no symmetry has been broken in the system, the state is 
smoothly connected to the regular SF state and no phase 
transition occurs in the usual sense. 

Diffraction Minima. Scattering light at 90°, one 
explicitly breaks the translational symmetry, which gets 
imprinted on the interaction of modes. In this case, we 
have two components in the diffraction minima in diago¬ 
nal density scattering reduces to ? 7 d,± = EJ^NgAp— 
with Ap = (p+ — P-)I2, where p± correspond to the 
mean atom numbers of each of the two light-induced 
modes. Thus, the effective chemical potential is p± = 
p ± gesJ^NsAp - gesJo and U± = U + 2gfeff J|,. In the 
absence of bond scattering {Jjj = (—I)^Jd Jij = 0 ), 
the system can support supersolid (SS), density wave 
(DW) insulators, MI and SF (Fig. [^). The system 
supports DW insulators present at half-integer fillings 
surrounded by a SS shell while Mi’s exist at commen¬ 
surate filling. There exists a critical light-matter cou¬ 
pling (|(/eff lA^s > U /2), where MI and SF are entirely sup¬ 
pressed, which is a quantum critical point. The system 
only supports supersolid (SS) and density waves (DW) 
due to incommensurability (Fig. [^). It forms giant 
density waves with maximal amplitude constrained by 
the mean atom number per site and is described by a 





































Devil’s staircase [IH]. Interestingly a DW insulator phase, 
checkerboard lattice, will appear above the giant DW 
transition at p = 2 , consistent with m- 

For off-diagonal coupling, when bond ordering occurs 
{Jij = {—lyjB Ji,i = 0 ) the situation is more subtle 
and this requires four components, 

t(pq = ^0 ~ ffeffPBjVS,; ( 24 ) 

9 ' = 1 

(25) 

where the component g -|- 4 is the same as q. This is the 
origin of dimer states discussed in the main text. 

The full entangled light-matter state. The full 
light-matter state can be found using an alternative pro¬ 
cedure to PEO] using rotations over the Hilbert space 
via canonical transformations [55]. We obtain 

(26) 

‘Pq 

where the subscript “a” (“ 6 ”) corresponds to the light 
(matter) part. The coherent state components for the 
light are: 


Wif^ja — D[a^py)\Q)a (27) 

with D{a) = exp(Q;d^ — a*a), the displacement oper¬ 
ator. The coherent state amplitudes are = cF^^, 
where c = g 2 /{^c F i^) ■ These amplitudes depend on 
the projection F^yipq)i, = F|(pq){,, that corresponds to 
the steady state component \^Pq)b of the effective matter 
Hamiltonian "Hgfj. The phase factor is due to the time 
evolution of the effective matter Hamiltonian: r^(t) = 
exp(-i'Hggt), with projection r^J(pq)t, = f^|(p,)t,. This 
is accurate in the limits |k/Ac| 1 , |to/Ac| ^ 1 , 
\JbEr/^ c\ ^ 1 - The ground state of the effective 
Hamiltonian is 14')^ = "Y^q^ \Fq)b^ where the number of 
components related to the light induced modes is f{R). 

For the multimode DW’s (i? > 2) originating from 
density scattering, the number of components f{R) > 
2 depend on the sub-lattice structure generated. The 
coherent state amplitudes are and depend 

on the projections Nq,yipq)b = D\ipq)i,. 

Similarly, for dimer states in mean-field approximaton, 
with 

BiPg = RlsJB[(l>qF4‘q+2—{<t>q+l+4>q+d)cOs{/S.(j))]/2, (28) 

while (j)q = = arg(V: 9 + 2 ) - a.ig{i{jq) = 

arg('!/:q+ 2 ) — arg('!/)g+i), and f{R) = 4. The components q 
and <7-1-4 are the same and ij^q = (bq) is the order param¬ 
eter. For SS dimers <j)q ^ 4>q+2 and (pq+i 7 ^ 4>q+^^ while 
for SF dimers (fiq = 4'q+2 and (pq+i = (l)q+ 3 - 


We see that both density and bond multimode self- 
organized states are represented by the superpositions 
of several macroscopic components. Using the quantum 
measurement of light, one can project the full states to 
one or several degenerate components, inducing novel 
states and dynamics |32|- 


Numerical methods. Numerical simulations of the 
effective Hamiltonian have been carried out by doing spa¬ 
tial mode decoupling theory truncated in Hilbert space 
using the Gutzwiller ansatz with self-consistent condi¬ 
tions with light-induced interaction given by Eqs. 0 
and (14). The problem is solved as a multimode con¬ 
strained nonlinear optimisation problem in 2R(f + 1) 
dimensions and 3i? constraints with R being the num¬ 
ber of light-induced modes and / is the filling factor of 
the mode occupation. For our parameter range it is 
sufficient to consider / = 5. Exact diagonalization cal¬ 
culations are done in the effective Hamiltonian Eq. (4) 
in the main text in ID for small lattices in the absence 
of tunneling. We calculate the largest eigenvalue of the 
reduced density matrix A scaled by the atom number N, 
X = X/N > 0. The condensate fraction estimated by 
A is an upper bound to the thermodynamic limit value 
of the SF fraction in the system [55| [5S] . Checkerboard 
insulators have A = 0 and SS has A 7 ^ 0. 
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